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ABSTRACT 



We present a model-independent technique for calculating the time of mid-transits. This technique, named "barycenter method" , 
uses the light-curve's symmetry to determine the transit timing by calculating the transit light-curve barycenter. Unlike the 
other methods of calculating mid-transit timing, this technique does not depend on the parameters of the system and central 
star. We demonstrate the capabilities of the barycenter method by applying this technique to some known transiting systems 
including several Kepler confirmed planets. Results indicate that for complete and symmetric transit light curves, the barycenter 
method achieves the same precision as other techniques, but with fewer assumptions and much faster. Among the transiting 
systems studied with the barycenter method, we focus in particular on LHS 6343C, a brown dwarf that transits a member of 
an M+M binary system, LHS 6343 AB. We present the results of our analysis, which can be used to set an upper limit on the 
period and mass of a possible second small perturber. 

Key words, planetary systems- stars, methods: data analysis 



1. Introduction 

The success of the transit-timing variation (TTV) method 
in characterizing p lanets around the stars Kepler 9 



X: 



m cn aractenzmg planets around tne stars Kepler y 
( Holman et al. 2010l ) and Kepler 11 ( Lissauer et al.ll2011 ) 
and i n detecting a planet around star Kepler 19 



amplitude of a TTV signal varies with the mass and dis- 
tance of the perturbing body. The latter has been used 
in several null detections to place an upper limit on the 
mass and or b ital parameters of a hypothetical per t urber 
feeaiJ 120091: ICsizmadia et al.l l2010l: [Adams et all |2010L 
201ll:l 



( Ballard et al.l 120111 ) strongly suggests that TTV method 
has come of age and is now among the main mecha- 
nisms for detecting extrasolar planets. This method which 
is based on modeling the variations that appear in the 
times of the transits of a planet due to the perturba- 
tions of other objects, has been shown by many au- 
thors to be capable of detecting small Earth-sized planets, 
moons of giant planets, and stellar c ompanions around 
varie t y of stars (|Miralda-Escuddl2002l; iHolman fc Murray! 
20051: lAgol et all 120051; iKippind 120091: iMontaltol 120101 : 



201ll : lMacieiewski et al.ll2010l l2011aMbT h 



The TTV amplitude is strongly amplified when 
the transiting and perturbing planets are i n a mean- 



motio n r esonance. For i nstan ce, as shown by lAgol et al 



(2005h: ISteffen et all 




lAgol fc SteffeJ (l2QQ7h: 



lHaghighipour fc Kirste 



Sch warz et al.ll201ll : lHaghigh ipour fc Kirstdl201ll ). 



Because, the interaction between the transiting planet 
and the perturbing body(ies) is gravitational if there are 
no magnetic fields, the amplitude of the TTV strongly 
depends on the masses of t hese objects and their or - 
bital architecture . As s hown by Holman & Murrav ( 20051 ), 
Agol et all (|2QQ5h . and lHaghighipour fc Kirstd (|201l| ). the 



lHaghighipour et al 
(|2Q1 lh a planet as small as Earth can produce large 
and detectable TTVs on a transiting Jupiter-like body 
in or near a resonance. This characteristic of resonant 
transiting systems makes the TTV method a powerful 
technique for detecting low-mass planets. 

The fact that different orbital configurations of the 
transiting and perturbing bodies can produce similar 
TTVs has made the inference of the mass and orbital 
elements of the perturber from the measurements of 
the transiting planet's TTVs a very complicated task. 
Several attempts have been made to overcome these 
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difficulties (jNesvornv fc Morbidellil 120081: iNesvorp |2QQgt 



Nesvornv fc Beaugel \20ld iMeschiari fc LaughliiJ l2Q10h . 
However, the complications still exist, particularly when 
the system is in or near a resonance . As shown by 
Garcfa-Melendo fc Lopez-Mor ales too continuous ob- 
servations by Kepler and C0R0T are expected to resolve 
some of these difficulties. 

Determining variations in transit timing requires pre- 
cise measurements of the times of mid- transits. To com- 
pute a mid-transit time, it is necessary to develop a theo- 
retical light-curve that best models the observational mea- 
surements of the intensity of the light of a star. When 
studying transiting planets, m any authors use the ana lyt- 
ical methodology developed by Man del fc Ago 1 (120021) for 
this purpose. In this method, the light-curve of a star is 
calculated using an analytical formula that contains sev- 
eral parameters such as the coefficients of the star's limb 
darkening, the ratio of the radius of the planet to that of 
the star, the semimajor axis of the planet (or its orbital 
period), and the planet's orbital inclination. To measure 
the individual mid-transit times, it is customary to hold 
all parameters (except mid-transit time) constant during 
the fitting procedure. As a result, the measurement of the 
time of each mid-transit will be vulnerable to systematic 
errors. In other words, any modification to the values of 
any of the above-mentioned parameters (which may be 
obtained when observing the system for longer times) will 
change the fitted light-curve and result in different values 
of the times of mid-transits. Subsequently, the values of 
the TTVs obtained in these systems will also change. 

We u sed a mo del- indepen dent methodology, first intro- 
duced bv lSzabo et al. ( 2006f ), to calculate the time of mid- 
transit. We call this technique the barycenter method be- 
cause it calculates the mid-transit times by using the defi- 
nition of the transit light-curve barycenter and its symme- 
try. We describe this methodology in section 2 and present 
examples of its application to some of the already known 
transiting systems in section 3. In section 4, we apply this 
technique to the system of LHS 6343 and explain its im- 
plications for the transit timing of the system. We analyze 
the derived O-C diagram of LHS 6343 in section 5, and 
in section 6 we conclude this study by summarizing our 
analysis and reviewing the results. 

2. Barycenter method 

As mentioned earlier, to determine the variations in the 
transit timing of a planet, a precise calculation of the times 
of its mid-transits is required. The mid-transit times are 
determined by fitting an analytically obtained light-curve 
to the observational data, and calculating the time of the 
mid-point of each individual transit on the latter curve. 
When the transiting body is planetary, the synthetic light- 
cur veJs_jis^iany_OToduce^ using the algorithm developed 
bv iMandel fc Agoll (|2QQ2h . In the majority of cases, the 
times of mid-transits are calculated by keeping all other 
parameters (e.g. stellar radius, planet radius, orbital pe- 
riod, and two limb darkening coefficients) constant during 



the fitting procedure. This is particularly important when 
the number of points inside a transit is small (e.g., 6-9 
points). In such cases, fitting the observed data can lead 
to imprecise results. 

Another technique for cal culating times of mid -transits 



is the "trapezoid method" ( Alonso et al.l l2QQ9h . In this 



method, a trapezoid function is fitted to the observational 
data and the best light-curve is determined by varying the 
depth, duration, and shape of the trapezoid. The time of 
mid-transit is then calculated by identifying the mid-point 
of each transit on the best-fit trapezoidal curve. 

In systems where the transiting/eclipsing body is a 
stellar companion, the time of each mid-transit /eclipse 
is calculated using the methodology developed by 
iKwee fc van Woerden (ll956l) . This method has been used 
by iDeeg et al.l ([2000L 120081 ) to calculate eclipse timing 
variations of eclipsing binaries caused by a circumbinary 
planet, and is based on the assumption that in an unper- 
turbed system, the light-curve of the transited/eclipsed 
star is symmetric. In this method, the mid-point of an 
eclipsing light-curve is determined by folding the light- 
curve around one point of the transit, and calculating the 
differences between the points on the two parts of the 
folded light-curve. The point where these differences be- 
come minimum corresponds to the point of mid-transit. 

Th e method developed by Kwee fc van Woerden! 
( 19561) has the adva ntage that unlike the method of 



Mandel fc Agoll (l2Q02h . it does not depend on the param- 



eters of the central star. However, for the measurements 
of the mid-transit times to be accurate, this method re- 
quires very many points, which are obtained through the 
interpolation of points from the results of observation. As 
a result, in transiting systems with few data points (e.g., 
the transiting systems identified in the long cadence of 
QO to Q2 data sets from Kepler), using this method is not 
practical. 

In this section, we explain a methodology that em- 
ploys similar idea as the method bv lKwee fc van Woerden! 
( 1956 ) (i.e., using light-curve's symmetry) and as such 
is independent of the system's stellar parameters. This 
method was first presented by ISzabo et al. and 
later used bv I Simon et al. ( 2007 ) and Kippind ~( 2011 ) to 



study the possibility of the detection of exomoons. It cal- 
culates the exact moment of mid-transit using the def- 
inition of the transit light-curve barycenter. We call this 
me thodology the "barycenter me thod" . Unlike the method 
by Kwee fc van Woerden! ([l956h , the barycenter method 
can be applied to transit planetary systems with few data 
points. 

To define the transit light-curve barycenter, we use a 
normalized graph of the flux of the central star. As shown 
in Figure 1, the flux of the star outside the transit is de- 
trended and normalized to 1. For a point i with a flux fa 
inside the transit, the corresponding value of the light-loss 
of the system is equal to 1 — fa. Similar to the barycen- 
ter point of a number of massive objects, we now define a 
barycenter for the points on the graph of the normalized 
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Fig. 1. Schematic view of the barycenter method. 
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Fig. 2. Comparison between the results of mid-transit of 
HAT-P-7b first transit, obtained by the Mandel & Agol 
flux. In this definition, the time of mid-transit will then method (black circle), the trapezoid method (red square) 



be given by 



and the barycenter method (blue triangle). 



£ - ft) 

1=1 

i=l 



(i) 



In equation (1), n represents the number (rank) of the 
data points in the observation, and t{ is Julian Day ( JD) of 
the observation point i. To obtain more precise results, we 
only consider the points that are inside the transit light- 
curve. A point is inside the transit if its light-loss (1 — fi) 
is higher than the standard deviation of flux outside the 
transit. 



3. Application of the barycenter method to 
known transiting systems 

3.1. HAT-P-7b 

To test the capability of the barycenter method and the 
validity of its results, we used this technique to calculate 
the times of mid-transits in several known transiting plan- 
etary systems. In addition to the timing, transits may also 
show variations in their durations and depths. However, 
our focus is only on the variations in the times of mid- 
transits. Our first case was the transiting planet HAT-P- 
7b (Kepler-2b). HAT-P-7 was observed in short and long 
cadences as a calibration target for Kepler. In the short 
cadence mode, the light-curve of HAT-P-7 consisted of ap- 
proximately 355 points in each transit. We used the results 
of the observations as reported in Q0 data set, and cal- 
culated the time of mid-transit for the first transit of this 
planet. Figure 2 shows the results (the first point from the 
left). The error bar on each point was determ ined using 
the bootstrap technique ( Wall fc Jenkinsll2003f ). As shown 
here, the time of mid-transit obtained from the barycenter 
method is consistent with thos e obtained from t he tra pe- 
zoid method and the model bv lMandel fc Agoll (120021 ). 



To evaluate the sensitivity of each of these techniques 
to the number of points in a transit, we reduced the num- 
ber of points in the light-curve by regular sampling, and 
calculated the time of mid-transit using all three methods. 
Results are shown in Figure 2. As expected, the sizes of 
the error bars indicating the uncertainties at each point in- 
crease for fewer data points. However, as Figure 2 shows, 
the times of mid-transits obtained by all three methods 
are close and agree with one another. 

3.2. Kepler-lb to Kepler-9c 

We also applied the barycenter method to the confirmed 
planets of the Kepler- 1 to Kepler-9 systems. Table 1 shows 
the results and their corresponding uncertainties. The un- 
certainties were calculated using the equation 



dT 



dfi 



(2) 



Table 1 also shows the values of the mid-transit t imes 
of these planets a s reported by iHolman et all (l201dh and 
Ford et al.l ([201 ll ) using the model of iMandel fc Agol 
(|2QQ2[ ). As shown here, the results obtained from the 
barycenter method agree very well with the previously 
reported values. 

Figures 3 and 4 show the differences between the 
values of mid-transit times ob tained by the baryc enter 
method and thos e reported by IHolman et al.l ( 2010f ) and 
Ford et aTl ( 2011 ). The error bar at each point was cal- 
culated by taking the quadratic sum of the uncertain- 
ties shown in Table 1. Table 2 lists standard deviations 
of these differences and their average error bars. As can 
be seen from Figure 3 and Table 2, Kepler-4b shows 
large error bars compared to those of other planets be- 
cause of its shallow transits. This figure also shows that 
the standard deviation of Kepler-3b, as listed in Table 
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Fig. 3. Differences between the values of mid-transit tim- 
ing obtained by the barycenter method a nd the Mandel 
fc Agol method for Ke pler- lb to Kepler-4b ( Holman et al] 
2Q10HFord etaHl201ll ). 



2, is larger than its average error, which can be at- 
tributed to the non-symmetric shapes of the first and 
sixth transits of this planet (see Fig. 5). These short-lived 
anomaly flux variations can be explained by different 
mechanisms such as the presence of a ctive regions (dark 

l2QQ9l 
2Qlll 



Silva-Vali 
201lh 



Jied a et a 
id I2QQ8I : 



spots) or a second transi t ing body (IRabus et al 
Sanchis-Oieda et all 120111 Sanchis-Oieda fc Win: 



Nutzman et al 



2011b iDeming et al 



Since in both the ba rycenter method and the method 
of lMandel fc Agoll ( 2002f ) it is assumed that the light-curve 
is symmetric, these methods are sensitive to missing points 
in the observation of a transit. This can be seen from 
Table 2 for Kepler-6b and Kepler-9b. The standard de- 
viations of these two planets are larger than their aver- 
age errors, which could have been caused by a missing 
point in the observation of the first transit of Kepler-6b 
and fifth transit of Kepler-9b (Fig. 6). To illustrate this 
effect, we made an artificial light-curve for a transiting 
planet an d calculated the t ime of its mid-transit using the 
model of Ma ndel fc Ago 3 (|2002l ). We then removed one 
point from the light-curve and calculated the time of mid- 
tr ansit using both t he bar ycenter method and the method 
of Mandel fc Ago 1 (l2002tl. As shown in Figure 7, and in 
agreement with Csizmadia et al. the results ob- 

tained by both methods show large deviations when the 
missing point was in ingress or egress (deviation ~ 500 
seconds). On the other hand, both methods become less 
sensitive when the missing point is close to the bottom of 
the light curve. This experiment suggested that both the 
barycenter and the Mandel & Agol methods require the 
full coverage of observation data in transit, and a missing 
point in the observation data may cause a large offset in 
the results. 

We also examined the applicability of the barycenter 
method to the long integration time of Kepler 's long ca- 



Fig. 4. Same as in Fig. 3 for Kepler-5b to Kepler-9c. 



dence observatio ns (29.42 minute s ) (iKippingl 201ol ). Using 
the algorithm by iMandel fc Agoll ([2002f h we generated an 
artificial light-curve with bins of long integration times 
(we chose a point every 6 seconds and used the mean of 
270 of those points as the observed flux). We changed 
the beginning time of each binning and studied the varia- 
tions of mid-transit times as determined by the barycenter 
method. Results point to a deviation of no more than 4 
seconds for the mid-transit times. 

4. The LHS 6343 system and its transit timing 

LHS 6343 is a close, M+M binary system with a separation 
of - 20 AU. The primary of this binary, LHS 6343 A (KIC 
10002261, RA=19h 10m 14.33s, Dec= 46 57' 25.50"), has 
a mass of 0.37 M and the mass of the secondary, LHS 
6343 B, is approximately 0.30 M (see Table 3). The pub- 
licly available Q0 and Ql data sets from Kepler revealed 
four deep transits in the light-curve of th is system. By 
analyzing these data, Ijohnson et al.l ( 201ll ) showed that 
these transits are produced by a third object, LHS 6343 
C, which orbits LHS 6343 A every 12.71 days. As deter- 
mined by these authors, LHS 6343 C is a brown dwarf 
with a mass of ~ 63Mj and is located at a distance 0.08 
AU from LHS 6343 A. 

In preparation for applying the barycenter method to 
the light-curve of LHS 6343, we analyzed each transit of 
this system separately. Our initial analysis of the light- 
curve of LHS 6343 at the time of the release of Q0 and Ql 
data sets pointed to a non-symmetric transit (transit num- 
ber 3) among the initial four transits of this system. We 
recall that the barycenter method is based on the symme- 
try of the shape of a transit. The release of the Q2 data 
set provided us with seven more transits, of which our 
analysis identified transits number 5 and 8 as asymmet- 
ric. Figure 8 shows these non-symmetric transits. We note 
that these anomalies may be caused by star spots. To bet- 
ter portray the anomalies in the shapes of these transits, 
we used Mandel & Agol's methodology and obtained the 
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Table 1. Transit timing of Kepler- lb - Kepler-9c, mea- 
sured by both methods, the first column is the number 
of transit, the second column lists mid-transits found by 
the barycenter method(M JD-24549 00), and the th i rd co l- 
umn lists mid-tra nsits reported in Holman et al. ( 2010f ); 
Ford et al.1 ([201 lh . 



Table 2. Comparison of standard deviation of the dif- 
ference between two methods and the average error bar 
obtained by quadratic sum. 



Number 


Tb-Barycenter method(M JD-2454900) 


T (M JD-2454900) 


1.1 


65. 645006±0. 000054 


65. 645036±0. 000073 


1.2 


68. 115621±0. 000051 


68. 115649±0. 000051 


1.3 


70. 586201±0. 000051 


70. 586262±0. 000079 


1.4 


73. 056852±0. 000055 


73. 056875±0. 000186 


1.5 


75. 527497±0. 000051 


75. 527488±0. 000098 


1.6 


77.998120±0. 000054 


77.998101±0. 000103 


1.7 


80. 468596±0. 000052 


80. 468714±0. 000071 


1.8 


82. 939301±0. 000051 


82. 939327±0. 000080 


1.9 


85. 409937±0. 000050 


85. 409940±0. 000074 


1.10 


87.880533±0. 000051 


87.880553±0. 000054 


2.1 


65. 381504±0. 000085 


65. 381496±0. 000185 


2.2 


67.586213±0. 000087 


67.586232±0. 000238 


2.3 


69. 791026±0. 000090 


69. 790967±0. 000178 


2.4 


71. 995973±0. 000088 


71. 995703±0. 000178 


2.5 


74. 200488±0. 000113 


74. 200438±0. 000181 


2.7 


78. 610067±0. 000089 


78. 609909±0. 000210 


2.8 


80. 814614±0. 000086 


80. 814644±0. 000172 


2.9 


83. 019473±0. 000086 


83. 019380±0. 000210 


2.11 


87. 428970±0. 000086 


87. 428851±0. 000199 


3.1 


67.591623±0. 000181 


67.587958±0. 000558 


3.2 


72. 476370±0. 000096 


72. 475763±0. 000343 


3.3 


77.363887±0. 000111 


77.363568±0. 000564 


3.4 


82. 250995±0. 000114 


82. 251373±0. 000555 


3.5 


87. 139357±0. 000103 


87. 139178±0. 000515 


3.6 


92. 026855±0. 000108 


92. 026983±0. 000491 


3.7 


96. 914806±0. 000094 


96. 914788±0. 000482 


3.9 


106. 690806±0. 000096 


106. 690398±0. 000368 


3.10 


111. 577701±0. 000096 


111. 578203±0. 000488 


3.11 


121. 353559±0. 000097 


121. 353813±0. 000448 


4.1 


104. 816415±0. 002636 


104. 815900±0. 002940 


4.2 


108. 029416±0. 002417 


108. 029584±0. 002460 


4.3 


111. 244774±0. 002342 


111. 243268±0. 003001 


4.5 


117. 672537±0. 002587 


117. 670635±0. 002792 


4.6 


120. 884243±0. 002438 


120. 884319±0. 002080 


4.7 


124. 098430±0. 002747 


124. 098003±0. 002296 


4.8 


127. 309430±0. 002302 


127. 311687±0. 002627 


4.9 


130. 527792±0. 002522 


130. 525370±0. 003034 


4.11 


136. 953717±0. 001792 


136. 952738±0. 002374 


4.12 


140. 167899±0. 001699 


140. 166422±0. 002767 


5.1 


66. 545846±0. 000776 


66. 546068±0. 000494 


5.2 


70. 094334±0. 000582 


70. 094536±0. 000520 


5.3 


73. 643932±0. 000517 


73. 643005±0. 000564 


5.4 


77. 191271±0. 000503 


77. 191473±0. 000434 


5.5 


80. 739494±0. 000514 


80. 739942±0. 000579 


5.6 


84. 288634±0. 000449 


84. 288410±0. 000543 


5.7 


87.836005±0. 000510 


87.836878±0. 000526 


5.8 


91. 385478±0. 000587 


91. 385347±0. 000495 


5.9 


94. 933849±0. 000492 


94. 933815±0. 000447 


6.1 


67. 430812±0. 000408 


67. 424550±0. 000441 


6.2 


70. 659594±0. 000359 


70. 659250±0. 000499 


6.3 


73. 894048±0. 000406 


73. 893951±0. 000421 


6.4 


77.128802±0. 000352 


77.128650±0. 000468 


6.5 


80. 363963±0. 000403 


80. 363351±0. 000453 


6.6 


83. 597843±0. 000355 


83. 598052±0. 000443 


6.7 


86. 832717±0. 000359 


86. 832752±0. 000568 


6.8 


90. 067911±0. 000355 


90. 067452±0. 000489 


6.9 


93. 302393±0. 000410 


93. 302153±0. 000495 


6.10 


96. 536655±0. 000356 


96. 536853±0. 000456 


7.1 


67.276389±0. 000429 


67.276027±0. 000374 


7.2 


72. 161758±0. 000406 


72. 161517±0. 000459 


7.3 


77.047387±0. 000407 


77.047008±0. 000382 


7.4 


81. 932553±0. 000332 


81. 932498±0. 000316 


7.5 


86. 817611±0. 000410 


86. 817988±0. 000455 


7.6 


91. 703986±0. 000412 


91. 703478±0. 000538 


7.7 


96. 588905±0. 000367 


96. 588969±0. 002094 


8.1 


64. 685860±0. 000691 


64. 686046±0. 000472 


8.2 


68. 208387±0. 000598 


68. 208545±0. 000345 


8.3 


71. 731198±0. 000581 


71. 731044±0. 000385 


8.4 


75. 253391±0. 000578 


75. 253544±0. 000447 


8.5 


78. 776414±0. 000555 


78. 776043±0. 000628 


8.7 


85. 821393±0. 000695 


85. 821041±0. 000602 


8.8 


89. 343196±0. 000618 


89. 343540±0. 000630 


8.9 


92. 866237±0. 000492 


92. 866039±0. 000644 


8.10 


96. 388414±0. 000731 


96. 388538±0. 000712 


9b. 1 


77.2484±0. 00086 


77.24875±0. 00087 


9b. 2 


96. 48276±0. 00099 


96. 4824±0. 00092 


9b. 3 


134. 95455±0. 00122 


134. 95437±0. 00077 


9b. 4 


154. 18997±0. 00094 


154. 19058±0. 00077 


9b. 5 


173. 44199±0. 00116 


173. 43412±0. 00107 


9b. 6 


211. 92629±0. 00082 


211. 92589±0. 00074 


9b. 7 


231. 17112±0. 00082 


231. 17167±0. 00071 


9b. 8 


250. 42960±0. 00074 


250. 42951±0. 00071 


9b. 9 


269. 68043±0. 00071 


269. 68103±0. 00068 


9c. 1 


69. 30489±0. 00091 


69. 30577±0. 00127 


9c. 2 


108. 33019±0. 00133 


108. 33086±0. 00111 


9c. 3 


147. 33572±0. 00107 


147. 3356±0. 00105 


9c. 4 


186. 31434±0. 00120 


186. 31251±0. 00107 


9c. 5 


225. 26218±0. 00103 


225. 26284±0. 00096 


9c. 6 


264. 18192±0. 00099 


264. 18168±0. 00100 



Planet 


Standard deviation (sec) 


Average error (sec) 


Kepler- lb 


3.3 


8.9 


Kepler-2b 


8.2 


18.6 


Kepler-3b 


76.0 


42.7 


Kepler-4b 


115.0 


306.3 


Kepler- 5b 


43.0 


65.0 


Kepler-6b 


168.2 


52.3 


Kepler- 7b 


26.5 


69.3 


Kepler-8b 


22.4 


71.3 


Kepler-9b 


233.2 


105.7 


Kepler-9c 


87.0 


132.7 



T able 3. Parameters o f the LHS 6343 system according 
to Uohnson et al.l ()201l[ ). 



Parameter 


Value 


Ma (solar mass) 


0.370 ± 0.009 


Mb (solar mass) 


0.30 ± 0.01 


Mc (Jupiter mass) 


62.7 ± 2.4 


a A B (AU) 


20.130 ± 0.605 


a A c (AU) 


0.0804 ± 0.0006 


Pc (days) 


12.71382 ± 0.00004 



best fit to all 11 transits of this system. Figure 9 shows the 
residuals of each single transit with respect to this best fit. 
As shown here, the residuals of the third, fifth, and eighth 
transits are larger than 1 — a (closer to 2 — a) because of 
their anomalies inside their transits. 

As mentioned for Kepler-9b, developing a model to ex- 
plain these anomalies would require many observational 
points inside each transit and will depend on several pa- 
rameters such as the size and latitude of starspots, their 
lifetimes, the rotational period of the star, and the orien- 
tation of the rotation axis of the star relative to the orbit 
of transiting brown dwarf. Although an interesting project 
on its own, developing such a model is beyond the scope 
of our study. Also, given that with the currently avail- 
able data, the number of points in each transit is limited 
to only 5 or 6, such a model may not even be entirely 
realistic. Therefore, because the basis of the barycenter 
method is on the symmetry of a transit, and also to re- 
strain false positive TTVs, we decided to exclude the most 
asymmetric transits (i.e., transits 3, 5, and 8) from our 
analysis. Table 4 lists the times of the mid-transits of the 
remaining eight transits of the system calculated using the 
barycenter method. To estimate the corresponding errors 
of each mid-transit tim e, we used the bootstrap method 
(|Wall & Jenkinsl l2003) and considered the standard devi- 
ation inside each transit as the initial uncertainty. Note 
that the standard deviation inside a transit may be larger 
than outside due to crossing over starspots. To check the 
validity of our error estimation, we also calculated the val- 
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A Sixth transit of Kepler-3b 
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Fig. 5. Asymmetries in the first and sixth transits of 
Kepler-3b (maybe due to starspots). 



First transit of Kepler-6b 
Fifth transit of Kepler-9b 



AAA 



Phase 

Fig. 6. Asymmetries in the first transit of Kepler-6b and 
fifth transit of Kepler-9b, caused by missing point of ob- 
servation. 

ue s of errors using equa tion (2) and the methodology used 
bv lDovle fc Deej (|2QQ4l ). Our calculations showed that the 
values of the errors obtained from all three methods have 
the same order of magnitude. 

To obtain the variations in the transit timing of the 
system, we applied a linear fit to the eight mid-transit 
times in Table 4. Results suggested a period of P = 
12.713815 days, corresponding to a semimajor axis of 
0.076-0.080 AU for the transiting body. These results 
closely agree with the results reported by Ijohnson et al 
( 201lh . 



Given that LHS 6343 is a binary system and the tran- 
siting object (LHS 6343 C) orbits the primary star, the 
three-body system of LHS 6343 AC-B forms a hierarchi- 
cal three-body system. We examined the stability of LHS 
6343 C in this system by numerically integrating its orbit. 
Results indicated that this object is stable for long times. 



400 



Barycenter method 
Mandel & Agol method 



1 23456789 

Num 

Fig. 7. Deviations of mid-transit timing from the known 
values as calculated by Mandel & Agol, and the barycenter 
methods for a synthetic light curve with one missing point. 
The x-axis presents the rank of missing points in light 
curve. 




Fig. 8. Strange non- symmetric shape of third (blue 
square), fifth (green triangle), and eighth (red circle) tran- 
sit of LHS 6343. 



We re fer the reader to a recent article bv lBorkovits et al 
( 201ll ) and references therein, where the authors have 
presented a detailed analysis of the dynamics and tran- 
sit/eclipse timing variations of hierarchical tripe systems. 

5. Analyzing the O-C diagram of LHS 6343 

The times of mid-transits obtained from the barycenter 
method show small deviations from their linear fit. Figure 
10 and the right column of Table 4 show these deviations 
and their corresponding uncertainties for each mid-transit 
time. In this section, we analyze these deviations from the 
linear fit of transit timings and discuss their implications 
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Table 4. Transit timing of LHS 6343, as measured by the barycenter method and O-C values in days were calculated 
according to the new linear ephemeris. 



Transit 


T (days) 


T (days) 


O-C (Sec) 


Number 


(The barycenter Method) 


(Calculated) 


(TTV) 


1 


54957.216473±0. 000133 


54957.216535 


-5.4±11.5 


2 


54969. 930434±0. 000154 


54969.930354 


6.9±13.3 


4 


54995. 358025±0. 000135 


54995.357992 


2.9±11.7 


6 


55020. 785698±0. 000120 


55020.785630 


5.9±10.4 


7 


55033. 499199±0. 000147 


55033.499449 


-21.6±12.7 


9 


55058. 927144±0. 000136 


55058.927087 


4.9±11.8 


10 


55071. 641056±0. 000173 


55071.640906 


13.0±15.0 


11 


55084. 354626±0. 000156 


55084.354725 


-8.6±13.5 



• Other transits 
■ 3rd transit 

A 5th transit 

• 8th transit 




0.00 0.02 

Phase 



Fig. 9. Residual of the best fit of the Mandel & Agol 
method to all transits of LHS 6343 (just inside transit). 




5 6 7 
Transit Number 



10 11 12 



Fig. 10. Diagram of the transit-timing variations of LHS 
6343. 



for the possible existence of a second smaller object around 
the primary LHS 6343 A. 

Because it is in a hierarchical tripe configuration, LHS 




At(scc) 



Fig. 11. Values of LTTs of LHS 6343 for different values 
of the binary semimajor axis. 



SxlO 4 



I0 4 



— — r 

0<c<0.S5 



3*10^ 



SxlO' 




3xl0 4 



10 4 



At (sec-) 



0,50<e<0,75 



3xl0 4 



2x10 4 



10 4 



i I i i ; 
0.75<e<0.90 



Arf sec) 



6343 C is continuously subject to the gravitational per- F i g . 12 . Values of LTTs of LHS 6343 for different values 
turbation of the secondary star. These perturbations af- Q f binary eccentricity, 
feet the orbit of this object and cause variations in the 
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Fig. 13. Graph of LTTs between 5 s and 6 s for different 
values of the binary semimajor axis and eccentricity. 



times of its transit (for a detailed analysis of TTVs in hi- 
erarchical tripe systems we refer the reader to Borkovits 
et al. 2011). Given that the semimajor axis of the binary 
(<~ 20 AU) and its projected separation (19-21 AU) are 
much larger than the semimajor axis of LHS 6343 C, it 
would be important to determine to what degree the vari- 
ations in the transit timing of this object have been caused 
by the binary's light-travel time (LTT) effect. To exam- 
ine this possibili ty, we used the methodology presented by 
iMontaftol ([20101 ) and calculated LTTs for different values 
of the semimajor axis and eccentricity of the binary. We 
changed the values of the projected separation of the bi- 
nary u sing the distribution given by iDuquennov fc Mayor 
(jl99lh . and performed 10000 LTT-calculations for ran- 
domly chosen values of the binary eccentricity between 
and 0.9. In all our simulations, we considered the system 
to be coplanar. We identified the systems for which the 
value of LTT was between 1 s and 6 s. Figures 11 and 12 
show the results for a timespan of three years (duration of 
Kepler's primary mission). As shown here, systems with 
LTTs between 4.5 s and 5.5 s constitute the majority of the 
cases (we recall that the minimum reported va lue of TTV 
that can be detected by Kepler is ~10 s see (IFord et al 



201lf )). Figure 13 shows the results of all our simulations 



for LTTs between 5 s and 6 s in more detail. As shown in 
this figure, the values of LTTs do not exceed 6 s, which im- 
plies that during the period spanned by the present public 
release of the Kepler observations data (~144 days), the 
contribution of the binary LTT to the variations in the 
transit timing of LHS 6343 C is negligible. 

The fact that the contribution of LTT to the O-C val- 
ues, as listed in Table 4, is negligibly small implies that 
in modeling these deviations from the linear fit of transit 



timings, one can safely ignore the effect of the secondary 
star. 

To examine whether these deviations could be caused 
by an additional body in the system, we considered the 
two-body system of LHS 6343 AC, and calculated the O- 
C values for different values of the mass, semimajor axis, 
and orbital eccentricity of a hypothetical perturber around 
the primary LHS 6343 A. To reduce the amount of calcu- 
lations (which could be large because of the large size of 
the parameter-space), we limited our study to only circu- 
lar and coplanar systems. Figure 14 shows the results for 
different values of the initial angular position of the hypo- 
thetical third body. A comparison between these results 
and the values of the O-C in the right column of Table 4 
suggests that a perturber with a mass ranging from 0.1 
to 1 Mj may be able to produce these values when in 
an orbit with a period ranging from ~ 3.5 to 8P (where 
P = 12.713815 days is the orbital period of LHS 6343 C) 
around LHS 6343 A. To determine an upper limit for the 
mass of the perturber, we calculated the O-C values for 
different values of the mass and semimajor axis of this ob- 
ject, and compared the results with the values of O-C as 
shown in Figure 10. Figure 15 shows the maximum val- 
ues of the mass of the perturber for which the value of 
X 2 between the O-C obtained from the model and those 
listed in Table 4 are lower than 3. As shown in this fig- 
ure, the mass of the perturber cannot be larger than one 
Jupiter-mass. 

6. Conclusion 

We presented a technique called the barycenter method for 
calculating the time of the mid-transit in transit-timing 
studies. This method is based on the symmetry of the 
light-curve, and has the advantage that is independent 
of the parameters of the system. In other words, unlike 
other techniques for calculating the mid-transit timing, 
the results obtained from the barycenter method will not 
change by changing the assumption on the parameters of 
the central star. However, the fact that this method re- 
quires symmetry in the light-curve implies that when the 
transit curve is not perfectly symmetric (i.e., when ob- 
servational points are missing, or because of starspots), 
large offsets may appear in the mid-transit timing mea- 
surements. The application of the barycenter method to 
several known transiting systems showed that the results 
obtained from this technique are comparable with those 
obtained from other methods. Our study indicates that for 
complete and symmetric transit light curves, the barycen- 
ter method achi e ves th e same precision as in the model of 
Mandel & Agoll (2002), but with fewer assumptions and 
much faster. 

We used the barycenter method to calculate the times 
of mid-transits of the M+M binary star LHS 6343. Our re- 
sults indicated that as suggested by Johnson et al (2011), 
the primary of this system is host to a smaller object with 
a period of ~ 12.7 days. A study of the variations in the 
transit timing of this body (LHS 6343 C) points to the 



10° 10 1 10 2 10 3 10 4 10 5 10 6 




2345672345678 



P2/P1 P2/P1 

Fig. 14. Values of TTVs in the system of LHS 6343AC 
due to a hypothetical perturber. The system is assumed 
to be circular and coplanar. Each panel shows TTVs for a 
different value of the angular phase of the perturber. The 
units on the color scale are in seconds. 
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4 4.5 5 5.5 6 6,5 7 7.5 8 

Fig. 15. Maximum mass of the perturber on circular orbit 
reproducing the TTVs as in Figure 10, as a function of its 
orbital period. 



possibility that a small object with a mass no larger than 
1 Mj may exist around LHS 6343 A, which can produce 
the O-C values lower than the upper values presented here. 
Whether such an object actually exists requires more tran- 
sit data and more observations of this system. 
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